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ABSTRACT 



The granulation pattern that we observe on the surface of the Sun is due 
to hot plasma from the interior rising to the photosphere where it cools down, 
and descends back into the interior at the edges of granules. This is the visible 
manifestation of convection taking place in the outer part of the solar convection 
zone. Because red giants have deeper convection zones and more extended atmo- 
spheres than the Sun, we cannot a priori assume that granulation in red giants 
is a scaled version of solar granulation. Until now, neither observations nor ID 
analytical convection models could put constraints on granulation in red giants. 
However, thanks to asteroseismology, this study can now be performed. The re- 
sulting parameters yield physical information about the granulation. We analyze 
~ 1000 red giants that have been observed by Kepler during 13 months. We fit 
the power spectra with Harvey- like profiles to retrieve the characteristics of the 
granulation (time scale Tgran and power Pgran)- We also introduce a new time 
scale, Teff, which takes into account that different slopes are used in the Harvey 
functions. We search for a correlation between these parameters and the global 
acoustic-mode parameter (the position of maximum power, fmax) as well as with 
stellar parameters (mass, radius, surface gravity (log g) and effective temperature 
(Teff)). We show that Tes oc ly'*^^ and Pgran oc which is consistent with 

the theoretical predictions. We find that the granulation time scales of stars that 
belong to the red clump have similar values while the time scales of stars in the 
red-giant branch are spread in a wider range. Finally, we show that realistic 3D 
simulations of the surface convection in stars, spanning the {T^s, log g')-range of 
our sample of red giants, match the Kepler observations well in terms of trends. 

Subject headings: red giants: general — methods: data analysis 
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Introduction 



Granulation was first observed by iHerschell ( llSOll ) on the Sun, and is widely known 
as the surface signature of convection, where bright cells of ascending hot gas are visible 
as the granules and the darker descending cool gas are the so-called "intergranular lanes" . 
These solar granules have typical sizes of about 1 Mm. The study of granulation is tightly 
related to the analysis of convection quantities, being the most important manifestation 
of convection at the surface of the Sun in terms of energy. Other phenomena related to 
convection are the acoustic (p-mode) oscillations that reveal the internal structure of the 
Sun. These oscillations are stochastically excited in the convective atmosphere. 

Red giants are cool, bright, and evolved stars. They have a surface gravity, log g, 
between 2 and 4 and an effect ive temperature, T^s, in the range ~ 4000 to ~ 6000 K, as 



defined by 



Ciardi et al. 



( I2OIII ). For these stars, the mass is r oughly in the 



4:Mq and log(L/L0) varies from 0.3 to 3 (IMiglio et al. 



2009 



Hekker et al. 



range 0.7 to 



2011af l. They 



are interesting not only because they provide constraints on di stance, age, and chemi cal 



Girardi &: Salaris 



200 ih but 



evolution of stars, galaxies, and the extragalactic medium (e.g. 
they also serve as laboratories for studying convection, granulation and p-mode oscillations. 
As in the Sun, the inefficient, super- adiabatic convection in the upper few pressure 
scale-heights is host to sonic, turbulent fiows that stochastically excite sound waves — these 
stochastic oscillations are also known a s solar-like oscillations . The unambiguous detection 



of non-radial oscillations in red giants ( iDe Ridder et al. 



20091 ) from photometric time series 



of several hundreds of red giants obtained with the CoRoT satellite (Baglin et al. 2006) 



was a significant step forward in red-giant seismology. These observations a . 



study p-mode glob al parameters and their scaling laws 



2010 



Mosser 



2010l ). to estimate their masses and radii (IKallinger et al. 



Hekker et al. 



2009 



retrieve evidence of sharp features in their internal structure (IMiglio et al. 



so allowed to 



Mosser et al. 



20 



Obf). a nd to even 



2OIOI) from, for 
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instance, the second ionization zone of helium. 



The launch of NASA's Kepler Mission (iBorucki et al. 



20101) in March 2009, took the 



next big step in asteroseismology. The Kepler field of view is 105 square degrees in the 
direction of Cygnus and Lyra and the filter of the photometer ranges from 4 300 to 8 800 A, 
with a broad peak at 5 900 A. In total about 150 000 stars are observed with high-precision 
photometry throughout the nominal lifetime (3.5 years) of the mission. While for ^17000 
red g i ants only a few mon ths of data have been released to the public domain, (jCiardi et al. 



2011 



Hekker et al.ll2011al ). around 1 500 red giants with time series longer th an a year are a t 



our disposal withiri the Kepler Asteroseismic Science Consortium (KASCl (IBedding et al. 



2010 



Hekker et al. 



2011b 



Huber et al. 



2010 



Kallinger et al. 



2010a 



Mosser et al 



2011ah . 



We thus have exquisite data for studying not only the acoustic modes and their dependence 
on stellar parameters, but we can also, for the first time, study properties of surface 
convection of a large sample of red giants. The observations and data processing that we 
have used are described in Section [2l 

When studying the p modes of stars, granulation is normally considered a "noise" term 
and referred to as "background". In the present paper, however, we use this background 
signal to investigate the surface manifestation of convection and the source of p-mode 
excitation. In Section [3] different methods used to extract characteristic granulation 
parameters from the power spectra are described and compared. In Section H] we investigate 
the correlations between granulation properties, stellar parameters (mass, radius, surface 
gravity and Tefj) and the frequency of maximum power, z/max- 

From a modeling point of view, convection can be treated either by ID analytical 



models, or by 2 or 3D hydrodynamical simulations. The t hree most emp 



models of convection ar e: the mixing- length theor y 



MLT formulations (e.g. 



Unno 



1967 : 



Gough 



1977a 



(MLT: 



Bohm-Vitense 



Dupret et al 



oyed 



analytical 



19581 ). non-local 



20061 ). and the Canuto 
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and Mazzitelli model (CM; ICanuto fc Mazzitellilll99ll . Il992[ ). All of these formulations of 
convection are employed in stellar-structure modeling, but have severe shortcomings in the 
surface lay ers, where most of the approximations they are built on, break down (for more 



details, see 



Trampedach 



20101) • In addition to these shortcomings, the free parameters of 
the analytical models suffer from large uncertainties. As a result, the analytical models have 
only limited predictive power. None of the analytical models deal with granulation and 
estimates of sizes, contrasts and flow speeds rely on a number of additional assumptions. 

Contrary to ID analytical models, realistic numerical simulations of convection are 
based on the quantum mechanics of the equation of state and opacities and the fundamental 
physics of radiative transfer and hydrodynamics. However, the simulations are also limited 
by the available computational power and hence employ various approximations to make 
them tractable. Most of the se approximations, s uch as the effect of a limited numerical 



resolution can be quantified (lAsplund et al. 



2000). 



Although 2D simulations can be run at much higher resolution, the topology of 
convection is fundamentally altered by the lack of a third dimension, and the properties 
of the 2D simulations cannot be directly related to observations of stars. Since the top 
boundary of convective envelopes occurs in the photosphere, a realistic treatment of 
radiative transfer is also necessary when the aim is to p roduce simulations that can be 



directly compared to observations ( iNordlund et al. 



2009 



, and references therein). The 



surface layers are, of course, the very layers we observe, which is the reason why we turn to 
more realistic, 3D convection simulations for guidance in the interpretation of the Kepler 
observations. 



Previously, 



Collet et al. 



(120071 ) used 3D red giant simulations for abundance analysis 



for [Fe/H] = to -3, u sing the same code and atomic physics as presented here 



(IStein fc Nordlund 



20001 ). They find metallicity effects on the granulation with increased 
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size of the granules for increased metallicity values. Some 30 rad ia tive-hydrodyriamics 



simulations are also done for subgiant stars f lRobinson et al 



20041 ) 



Ludwig et al 



fl2009[ ) 



compared their 3D radiative-hydrodynamics simulation carried out with CO^BOLD code 
with the CoRoT observations of the F-type dwarf, HD 49933. They found a significant 
over-estimation of the theoretical signal by a factor of two to three in total power. 



Stein fc Nordlundl (120001 ) code 



In this work, we use the simulations performed with the 
and compare them to observations of ~ 1000 red giants. The simulations are introduced in 
Section 5, where we also explain the synthetic Kepler data of these simulations and compare 
the results of the simulations with the observations. 



2. Observations 

2.1. Sample selection 

In the present investigation we used Kepler observations of red giants showing solar-like 
oscillations, which were pre-selected for asteroseismic or astrometric purposes before launch 
of the spacecraft. While the asteroseismic sample is composed of subsamples with different 
selection criteria, the astrometric sample, which consists of the majority of our total sample, 
has been selected coherently. These ~1 000 stars are selected to be 1) as distant as possible, 
to ensure small parallaxes and proper motions, 2) as bright as possible, but with very 
small chances of saturating the detector in any season, 3) located in fields as uncrowded as 
possible, and 4) uniformly distributed over the focal plane and not close to the edge of a 
CCD. This resulted in the following criteria: 

• Teff < 5400 K 

• log^f < 3.8 (c.g.s.) 
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11.0 mag < Kepler magnitude < 12.5 mag 



crowding > 0.95 (i.e. very low probability of having light from another source) 



Tpff and log q 



2010 



Brown et al. 



Tom the Kepler Input Catalog (KIC, iBatalha et al.ll2010l : iKoch et al 



20 111 ) have been used. 



Long-cadence (~29.4 minute sampling) data of the first ten days (QO), the initial 
one-month roll (Ql) and each of the following three-month rolls (0,2, 03, 0,4, and Q5) h ave 



been used to obtain long time series of more than one year of da ta (Jenkins et al. 



We have used data processed as described in 



Garcia et al. 



2010|). 



(I2OIII ) in which instrumental 



effects, such as satellite safe-mode events or attitude adjustments, have been removed. 
These effects can affect the high-frequency range when jumps happen in the light curve or 
the low-frequency range when a safe-mode event is observed. Long-term variations have 
been preserved, which could contain signal of stellar origin, e.g., granulation, but also 
from instrumental effects. The latter are filtered out with a procedure described below 
(Section 2.2). 

Among the 1283 stars for which we have Q012345 data, we have selected a sample 
of ~ 1000 red giants (900 to 1100 depending on the method used to analyze them). 



corresponding to the ones where different teams returned values in ag reement for t 



ne mean 



large separation, (Ai/), and the position of maximum power, Uraax- See 



Hekker et al. 



(12011bh 



and references therein for descriptions and comparisons of the different methods and results. 



2.2. Removing instrumental long term variation 



2OIII ). we need to 



From the long-term variations preserved in the data (1 Garcia et al. 
remove the instrumental signal without affecting the granulation signal. We assume that 
the characteristic granulation time scale Tgran is of the same order as the time scale of the 
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oscillations, as is the case for the Sun (see e.g. I Vazquez Ramio et al.ll2005l ). For red giants 
we therefore expect Tgran to be of the order of hours up to a few days for the largest stars 
of our sample. Note that these equivalent time scales differ by roughly a factor of vr in the 
frequency domain as an e-folding time of r in power equates to a position in the frequency 
domain of l/(7rr), while a wave of period T has a frequency of 1/T. 

We have investigated effects of different high-pass filters used to remove the signature 
of signals with periods longer than the expected Tgran for red giants, which we assume to be 
(partly) due to instrumental effects. We tested different filter shapes and decided to use 
a triangular smooth as it introduces less ripples due to sharp edges compared to boxcar 
filters. Additionally, we tested different filter widths of 10, 20, and 30 days. It showed 
that the wider filters result in increasingly different Tgran values obtained from the data, i.e. 
introducing a bias as well as additional scatter. From these results we concluded that a 
triangular smooth with a full width at half maximum of 10 days applied to an interpolated 
time series provides most robust results. Note that the filter with a width of 10 days 
limits our sensitivity of Tgran to periods < 10 days, which is sufficiently long to study the 
granulation in red giants under the aforementioned assumption that these are of the order 
of hours up to a few days. 



3. Granulation parameters 



HarvevI fllQSSf ). He 



Granulation parameters were first modeled in the Sun by 
approximated the autocovariance of the time evolution of the granulation by an exponential 
decay function with a characteristic time Tgran- This results in a power spectrum with a 
Lorentzian profile: 



1 + (27rz/Tj 



gran 



(1) 



grany 
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in which Pnii^) is the total power of the signal at frequency u, a is the characteristic 
amplitude of the granulation and a is a positive parameter characterizing the slope of the 
decay. We also define the amplitude of the granulation power, Pgran = 4cr^rgran- This 
approach has been successfully applied to the solar m easurements by several instruments in 



both velocity and intensity (e.g. 



Lefebvre et al. 



20081). 



Other signals such as mesogranulation, supergranulation and active regions can be 



modeled with similar 



mentioned in [Harvey 



uncti ons, and act at much longer time scales. We note here that it is 



19851 ) that the declining slopes of the functions may well be different 



from two. Di f ferent modified Harv e y fun c tions have in d eed b e en used by for instance 



Aigrain et al. 



2ooi 



Hekker et al. 



(120 loh . 



Huber et al. 



(12009f ) 



Kallinger et al. 



(12010ar i 



and in several methods used in this work, which are explained below in more detail. The 
use of values different from a = 2 were established empirically, but the physical meaning of 
the function with a=2 has not been deeply investigated previously. 

Recalling that the autocorrelation of a signal is nothing but the Fourier transform of its 
power spectrum, we can easily investigate the influence of a on the autocorrelation function 
(ACF) of the temporal signal of the granulation component. We note that, for a given r, an 
increasing value for a indicates an increase in the temporal correlation. For a = 2, the ACF 
is obviously an exponential function decaying over a time scale r. When a tends to zero, 
the power spectrum becomes flatter and flatter, i.e. tends to white noise. The ACF does 
then decrease faster and faster, and tends to a Dirac function, typical for the ACF of white 
noise. With increasing a, the ACF broadens and converges to the sine function for a — > oo. 

To be able to compare characteristic time scales r obtained with different a values or 
even with different expressions of the granulation spectrum, we define an effective time scale 
Teff as the e-folding time of the ACF. For a = 2, we thus recover r = t^s- In the following, 
we will study Tes instead of Tgran to correctly account for slopes that differ from two. 
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3.1. Methods 



Six teams have fitted the granulation background of the red giants and the values of 
the parameters have been com pared to check tlie vahd ity of the resuhs. For all the methods 



explained below, except CAN (iKallinger et al 



2010aJ), the fit is performed by not taking 



into account the frequency region in which the solar-like oscillations are present. CAN 
fits the p-mode bump together with the background. The methods are all summarized in 
Tables [Hand 12 

The OCT method fits the background of the power density spectrum (PDS) using one 
Harvey-like model: 



(2) 



The initial input parameters are [Pgran.o, Tgran,o, «o, Wq] = [P, 0.01, 2, noise], in which P is 
the maximum power of the binned power spectrum (bins vary with frequency and are 2Au 
wide), and noise is the minimum power of the binned power spectrum. All four parameters 
are left free in the non-linear least-squares fit. When 100 iterations are done, the model 
with the lowest is chosen. The standard dev iations of the fitted parameters are used as 



the uncertainties. More details can be found in 



Hekker et al. 



teoioh 



The method adopted by SYD (see 



Huber et al. 



background with two modified Harvey functions (IKaroff 



2009 



, for details') consists of fitting the 



20081 ) : 



PM 



i=l 



(3) 



with Oi the root mean square (rms) intensity in power density, ri corresponds to the 
granulation time scale and T2 might correspond to another feature with no physical 



13 



meaning so far. This method uses a Levenberg-Markward least-squares algorithm to fit the 
background and to estimate Uj and Tj. The initial guess value for the granulation time scale 
is scaled from the Sun using the relation Tgran oc l/z^max, where i^max is the frequency of 
maximum oscillation power and is determined before the background model is fitted. To 
stabilize the fit, only frequencies in the power spectrum lower than 2z/inax are used for fitting 
the granulation component. This fit is performed on a binned PDS, leading to a Gaussian 
distribution of the points in the power spectrum. Monte-Carlo simulations are used to 
estimate uncertainties on the fitted parameters. For each simulation, a synthetic spectrum is 
generated by drawing random values following a distribution with two degree of freedom 
(d.o.f.) and repeating the fitting procedure described above. The uncertainties are taken as 
the standard deviations of the distributions derived from typically 1000 simulations. 

The CAN method uses three Harvey-like functions to model the stellar background 
signal and a Gaussian to account for the additional power due to p-mode oscillations: 



with the condition ti < r2 < T3, where the indices indicate consecutive background 
components. The Gaussian p-mode component has amplitude Pg, central frequency 
z/jiiax and width Ug. This method uses a Bayesian Markov- Chain Monte-Carlo (MCMC) 
algorithm. 

The parameter T2 corresponds to the granulation time scale. For now, a physical explanation 
of the other time scales, t\ and T3, is lacking. The latter could be related to a different 
scale of granulation and is currently under investigation. The granulation frequency, 
1/T"gran = ^1^2-, was allowed to vary from to 10 times t'max- -Pg was allowed to vary from 
zero to 10 times the average power in the spectrum around the initial guess for i^max, and 
W was kept between 0.5 and 2 times the average power at the high frequency end of the 
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spectrum. The adopted likelihood function is based on a distribution with 2 d.o.f. 
and the uncertainties are determined from the marginalized posterior probability density 
distributions. 

The DLB method (D. Buzasi, private communication) fits a single component Harvey 
model plus white noise as indicated in Eq. |2]with a = 2 as a fixed parameter. This method 
uses a linear least squares fit over the frequency range [2/T, z/Nyquist], wherer T is the 
length of the time series. Initial guesses for the fit [W, Pgran, Tgran] are given by [Pmin, Pmax, 
3 X 10^] , where Pmin and Pmax are respectively the minimum and maximum power found in 
the region considered in the PDS. No uncertainties have been estimated with this method. 

The COR method (B. Mosser, private communication) focusses on the frequency range 
around the break of the Harvey functions, just below the identified oscillation power. The 
fit of the background is performed with a single Harvey-like function in the frequency range 
between z/max/30 and 2z/ma,c- The initial guesses of the fit are derived from the mean values 



of the background at u^ 



region, 



Mosser et al. 



(IMosser et al. 



2nilc 



aj) and at low frequency (in the plateau 



20101). The exponent of the Harvey component (Eq. 1) is a free 



parameter. 



With the A2Z metho d, the background is fitted with one Harvey-like model and a 



power law as described in 



Mathur et al 



(120101 ): 



P(iy) = W + Puiiy) + au-\ (5) 

The initial guess for W is derived from a general relation between the stellar magnitude 
obtained from the KIC and the noise in the time series. This relation was determined from 
a larger sample of solar-like stars observed in short cadence by Kepler. The initial guess of 
the exponent a is 2, while that of Tgran is set to 2.1 xlO*^ s. Assuming a distribution 
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with 2 d.o.f., the fit is performed using a maximum hkehhood estimator. The uncertainties 
are estimated by inverting the hessian matrix. 

For each method, we have computed Tcs converted from Xgran and Pgran as defined in 
Eq. (1). 

3.2. Comparison of the results 

Figure [1] shows a typical PDS of a red giant observed by Kepler. The results of each 
method listed above are shown and all of them provide the same qualitative results around 
the break of the granulation profile with small differences at very low frequency or at very 
high frequency. We note that one method can work well for one star while another method 
might work better for another star. 

Figures [2] and [3] show the results of r^s and Pgran for the methods described in Section 
3.1. For the methods where more than one Harvey-like function was fitted, the one with 
the smallest time scale has been interpreted as belonging to the granulation for the SYD 
method and the middle time scale was used for the CAN method. The trends in the results 
from different methods are consistent, although they are not always the same. Interestingly 
the correlation between t^s and Pgran as shown in Figure H] is more coherent for the different 
methods. This could indicate a degeneracy in the fitted parameters coming from the use 
of different methods with different free parameters, models, data, number of components 
as it can be seen in Figure [TJ However for a given method, the results and the trends are 
consistent. 

For the different methods, the average uncertainties on l/r^s are 3%, 9%, 9%, 1%, 
and 6% for A2Z, CAN, SYD, OCT, and COR respectively. Some of these uncertainties (in 
particular for A2Z, OCT and COR) seem rather small and might be underestimated so we 
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assume that the uncertainty should rather be around 10%. If we look at the dispersion 
between the methods for a given Uma.x, 1/Tefr varies by ~ 30% around the average value. 
The dispersion between the stars can be related to the different metallicities of the stars. 

One thing that is apparent is the existence of a "jump" in the Tgran results at 
i^max ~40 /xHz for two of the methods using a single power law, an initial guess for Tgran 
independent of Vmax and a as a free parameter. This second branch is most likely due to 
the fact that the fit converges to a second feature in the power spectrum. The absence of 
the second branch in the results of DLB (guess of Tgran independent of z/max but fixed a) 
indicates that the convergence is at least partly due to the fact that a is a free parameter. 
We also performed tests with different initial guesses for Tgran (again independent of i^max) 
and that reduces the second branch but the fit does not converge for high values of t'max- 
Additionally, to investigate further the number of components, we performed tests with one 
method, A2Z in this case, in which we fitted either one or two Harvey components with 
several initial parameters to investigate how these affect the results (see the appendix for 
more details). It appears that by using two Harvey-like functions, we obtain only one branch 
instead of two branches (see discussion in Section 6). For all of these reasons, we think 
that on the one hand, the minimization algorithm of the fit could go to a local minimum 
instead of a global minimum, which can lead to the second branch. On the other hand, this 
second branch could be related to the initial guesses, for which we converge to a background 
component with lower Tgran, which is likely to be something different from granulation. 
Prom the fact that the second branch also shows a strong correlation with i^max, we expect 
it to be of stellar origin and it could correspond for instance to a phenomenon similar to 
bright points or stellar spots whose appearance is modulated by the stellar rotation. 

As the results of the methods are qualitatively similar, we use one method for the 
representations of our results in the remainder of this work. 
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4. Granulation and other stellar parameters 



4.1. Granulation scaling relations and observations 



From some basic physical assumptions, it is possible to predict the b e 



lavior of 



granulation with respect t o stellar parameters. As argued by 



Huber et al 



(120091 ) and 



Kjeldsen fc Bedding] (j201l[ ). the convection cells travel a vertical distance that is proportional 



to the pressure scale height, Hp, at a speed that is approximately proportional to the sound 
speed, Cs- Since Hp oc Tes/g and Cg oc y/T^, the characteristic time scale of the granulation 
can be expressed as Tgran oc Hp/cg oc y/T^/g oc y/T^R^ / M oc L/{MT^^) oc 1/z^max, where 
g is surface gravity. Teg is effective temperature, R is stellar radius, and M is stellar mass. 
This result indicates that we should expect longer characteristic granulation time scales on 
larger stars, i.e. st ars with lower surface gra vity. Support for the argument is given by the 



3D simulations by 



Trampedach et al. 



(120111 ) carried out with the Istein fc Nordlundl (I2000l ) 



code. These simulations indicate that red giants have granulation cells that are roughly 



15 Hp wide. 



The proportinality of the horizontal size o 



H 



(e.g., 



Schwarzschild 



1975 



Antia et al 



a granule, d to the pressure scale height. 



19841 ) implies the following: d oc Hp oc T^s/g oc 



TesR'^/M oc VTeff/i/niax- Thus, the linear size of the granules is inversely proportional to 
the gravity, and since Tes and the mass of the star vary less than the radii of red giants, it 
indicates that the size of the granules should increase with the radius of the star. 

For the number of granules on the surface, N, we derive the following: oc R'^/d'^ oc 
M^/(Tjji?^) oc Mg/T^g. This indicates that the number of granulation cells decreases with 
increasing stellar radius or equivalently decreasing surface gravity. So, large stars have a 
relatively small number of large granulation cells on their surfaces compared to smaller 
stars with a relatively higher number of smaller cells. These results are in line with earlier 
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work by e.g. ISchwarzschildl (Il975[ ): lAntia et al.l (Il984j ). Evidence for large granulation 



features has been o 
review see 



2009 



j served on red supergiants such as Betelgeuse usin g inter f erometry (for a 



Monnier 



20 



Kiss et al. 



(120031 ) ■ for recent observations see 



Haubois et al 



1 ) and for their interpretation with 3D simulations, see 



(20091 



Ohnaka et al. 



Chiavassa et al. 



mm; 



( I2OIOI )). We can test this hypothesis here also for red giants. 



As already pointed out by the original Harvey models, the power of the granulation 
-Pgran is proportioual to cr^Tgran, CT being the rms intensity fluctuation. The proportionality 
with Tgran also meaus that Pgran oc d, i.e. fewer larger cells result in higher granulation 
power, and thus each cell has a larger influence on the total luminosity fluctuations of 
the star, and hence a larger influence on its variability. This has some similarity with 



(Mosser et al. 


2010; 


Huber et al. 


2010) 



In Figures [2] and |3] we fitted a power law for each parameter and for each method. The 
values of the slopes are listed in Table [3l By taking the results of all the methods together. 



while Pgran OC t'max , which is close to the re 



we find that r^s oc z/, '^■^^ 
(Tgran 0^ '^max) with the relations obtained by 



ation derived above 
Kieldsen fc Beddind fl2011l ) (Pg^an oc u'^^J 



Finally, it also agrees with analytical models of convection where a coefficient —1 is found 
for the Tgran-^max relation (R. Samadi private communication). Given that both p modes 
and granulation are driven by convection, this correlation is not surprising. 



Figure m shows the variation of Pgran with Teg. The fit of a power law (see Table [3]) 
gives a slope of 2.18 for all the methods together, which agrees with the scaling laws. 

It is interesting to study the variation of the ratio between Pgran and t^q, which is 
proportional to the variance of the intensity variations. This quantity varies as a function of 
log g. This is illustrated in Figure |5l where we used the granulation parameters obtained by 
A2Z as an example, but this tendency is qualitatively similar for the other methods. Stars 
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that have smaller log g, thus larger radii present a higher intensity contrast compared to 
stars with a higher log g. This anti-correlation betwee n Pp-ran / Tpff and log q seen m Figure O 



is mostly due to the fact that Pgran scales with ^ (IKjeldsen &: Bedding 



20111 ). which 



arises from the averaging of fluctuations over many (unresolved) granules, diminishing the 
power for large A^, i.e, having fewer granules leads to less averaging over the stellar surface 
and the effect on luminosity from each granule is higher. 



4.2. Granulation and global stellar parameters 

We investigate the dependence of the granulation parameters of the red-giant stars 
with stellar fundamental parameters R, M, T^ff, and lo q. The se parameters have been 



computed with the method described by 



Kallinger et al 



teOlOal ) using 13 months long time 



series for a sample of 1035 red giants. Figure |6] shows how Tefj computed with method CAN 
varies with R, logg, M, and Tgfj. The results from the CAN method are shown here because 
this method provides results of a large sample of stars. As predicted from scaling relations 
(see Section 4.1, Tgran VT^R"^ /M) for bigger stars, the granulation time scales are larger 
and we can fit a power law showing the correlation between Teg and R as oc r^g. This is 
also consistent with the predicted anti-correlation between Tgran and logg. Correlations of 
granulation parameters with mass and radius are not as tight as with q = GM/R"^, since the 



Gai et al. 



20111 ). The correlation with 



latter is the underlying independent variable (see also 
Teff seems much less tight than with logs'. However, the plot of Tcs vs. Tcs is dominated 
by stellar evolution effects and hence shows a similar behavior as seen in the H-R diagram 
(Figure [8]) including the spread in points intrinsically linked to the distribution of stars in 
the H-R diagram. From the fact that Pgran oc Tgran, we would expect similar correlations for 
the granulation power. These are indeed observed and shown in Figure [3 



We note a bump around 10 P©, log g = 2.5 and Tes = 4700 K in Figures |6] and [71 It 
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corresponds to the so-called red clump, corresponding to red giants that have gone through 
a Helium flash and are now in their He-core burning phase. In the plot showing Teg as a 
function of T^s, we also clearly see the red-giant branch. 



4.3. Granulation in the H-R diagram 



Having a large number of red giants that are at different stages of their evolution allows 
us to study the distribution of the granulation parameters in this part of the H-R diagram. 

The left panel of Figure |8] shows the distribution of l/rgfj in our sampl e of red giants. 



As prev iously, the values of T^s and L have been computed as described by 



Kallinger et al. 



( l2010al ). For stars with increasing luminosity, l/r^s decreases, thus is larger, which 
also correlates with larger granulation cells. There is a large concentration of points at 
logL/IvQ ~ 1.7. These stars belong to the red clump. 

The distribution of the granulation power, Pgran, is shown in the right panel of Figure El 



We note that low-luminosity stars ha ve low P^a.„ va 



we observe for the mode amplitudes (IMosser et al. 



ues, which i s once again similar to what 



2010 



2011aJ). This is not surprising as 



the p modes are excited in the convection zone so the same order of magnitude of energy 
would be involved in both granulation and p-r aode excitation. This is in agreement with 
the study of the height-to-background-ratio by iMosser et al. 



(1201131 ). 



Recently, 



Beck et al. 



(1201 If ) 



Bedding et al. 



(120 111 ), and 



Mosser et al. 



fl2011b[ ) detected 



for the first time mixed modes in the red giants. They showed that the period spacings 
between the mixed modes make it possible to distinguish between stars ascending the 
red-giant branch, i.e., H-shell burning stars and stars in the red clump, i.e.. He-core 
burning stars. Using this information, we investigated whether the stars show differences 
in granulation parameters related to their evolutionary phase. Due to the uncertainties 
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in the granulation parameters a firm conclusion can not yet been drawn, although some 
indications of different granulation parameters with evolutionary state are present. It is at 
least clear that the red-clump stars have granulation values in a specific range, and the 
red-giant branch stars have their values spread over a larger range. It is not surprising as 
the structure of red-clump stars is very similar, which would lead to similar granulation 
characteristics. 



5. 3D Simulations of Convective Red Giant Atmospheres 



To help us interpret the Kepler observations, we have turned to 3D hydrodynamic 
simulati ons of convection in stella r surface layers. These simulations, described in more 



detail in 



Trampedach et al. 



(l201l[ ). were constructed for direct comparison with observations 



and are therefore based on the state-of-the-art atomic physics. Th e thermodynamics is 



provided by tables of the so-called MHD equation of state (EOS) (IHummer &: Mihalas 



1988 



Mihalas et al. 



19881 ). which is based on explicit accounting for all excited states 



in all ions of the 16 most abundant elements, and a physical model of the (non-ideal) 



effect of interactions 



jetween particles . The opacities are based on the Marcs stellar 



atmosphere package (iGustafsson 



19731 ). but with several upd ates of 



an d line-opacities frona the Atlas stellar atmosphere package ( iKurucz 



he opa city sources 



1992a 



m 



b|), as detailed 



Trampedach et al.i (120111 ). The three main processes governing the surface layers (i.e.. 



what we can observe) of late type stars are hydrodynamics, thermodynamics and radiative 
transfer. These processes interact in very complicated ways, hence the need for realistic 
first-principles simulations. The surface layers (by definition) straddle the photospheric 
transition from the optically deep interior and the optically thin atmosphere, which means 
the radiative transfer cannot be treated in any known approximation. The full wavelength 
and ray-direction dependent problem has to be considered and the simulations are based on 
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the opacity binning formulation introduced by iNordlundl ( 119821 ) . The simulations are each 
performed on a grid of 150 x 150 x 82 points, spanning about 13 pressure scale- heights in 
the vertical direction (with 7 below the photosphere) and reaching up to an optical depth 
of log^o ''"Ross = —4.5. The horizontal extent primarily scales with gravity, with only a slight 
increase with Tes and is large enough to cover about 10 major granules. Each simulation 
was relaxed to a statistically steady state, with no systematic drifts in fluxes or mean 
structure over time. Surplus energy was extracted by artificially damping radial p modes 
during the relaxation. After this relaxation, production runs were carried out, covering at 
least 10 periods of the fundamental p mode. The resonant modes of the box (the bottom is 
a node) are excited and damped by the convection in the box and saturate at an amplitude 
given by the balance between the damping and driving. With these fairly short time-series 
ranging from 14h55 for the coolest subgiant to 5d 17h30 for the hottest giant, we see two 
to three radial modes and similar for the non-radial modes that has the box-width as 



ho rizontal wavelength . AH 3 7 simulations of the grid were perforn ied for solar meta. 



m 



li city a s 



Anders fc Grevessd (119891 ) with He and Fe adjusted to mimick iGrevesse fc Noeld (119931 ) 
and cover the zero-age main-sequence from T^s = 4 300 K to 6 900 K and up to giants of 
logg = 2.2 between T^s = 5 000K and 4200K. Our sample of 1035 Kepler red giants is 
spanned by seven of the 37 simulations of the grid. About 11% of the stars, however, fall 
outside the simulation grid and are therefore not included in the comparisons of Section 15.21 

Similar simulations that are not part of the grid, have bee n widely applie d to so lar 



Asplund et al. 



( 120001 ) 



and stellar observations. A study of solar Fel and II lines by 
showed remarkable agreement with the shape, asymmetry and wavelength-shift (without 
adjustable parameters to match these profiles to the observations). Note that ID solar 
atmosphere models cannot reproduce these observations. A similar analysis was carried 
out for Procyon with a similar level of success. This agreement with detailed line-shapes. 



leads to a recognition of the Ni I blend with the [O] I line, leading to agreement between 
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the oxygen abundances from [O] I and O I-lines ( iPrieto et al 



200 



variation of Na- and 0-lines in the Sun also matches observations (jPrieto et al 



)■ The center-to- limb 



20041 ) ■ Even 



earlier simulations could repr oduce the observed siz es, shapes and contrast of granules 



( iNordlund fc Stein 



199l|) and 



Rosenthal et al. 



(119991 ) showed how the stratification change 



from a realistic 3D convective atmosphere accounts for most of the so-called surface term 
in helioseismology. Together, these disparate observational tests probe a large range of 
depths in the atmosphere of the simulation, giving us confidence that they are a better 
representation of stellar atmospheres than are ID models. Hence, the 3D simulations have 
a larger predictive power. 



5.1. "Observing" the Simulations 

For each simulation, radiative transfer was performed for the full set of wavelengths 
in the opacity distribution functions (ODFs). The ODFs that were used are defined for 
1100 wavelength regions, with 12 points in each of these distribution functions. We have 
artificially assigned wavelengths to these 12 points within their respective wavelength 
regions, in a way that gives them the appropriate integration weight. The resulting spectra 
are saw-tooth shaped with 20 A wide bins in the optical. This is obviously not adequate for 
monochromatic studies, but this wavelength re-ordering of opacity on a 20 A-scale, has very 
little effect on broad-band colors. Figure |9] shows such an ODF spectrum together with the 
Kepler filter. Having computed complete spectra for each snapshot of each simulation, and 
for eight /x-angles (/i = cos^ of the position angle on the stellar disk) and four azimuthal 
0-angles, we then proceeded to average over and convolve with the transmission curve of 
the Kepler fAter. This results in Kepler intensities as function of /i-angle and time, I{fi,t). 

The subsequent transformation f rom specific intensities, to power spectra of observable 



fluxes, was performed as described by 



Trampedach et al 



(Il998l ) and later described in great 
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detail by iLudwid (|2006[ ). Following the derivations of iLudwia (|2006[ ) we write 



(6) 



1-3 



where I is the horizontal extent of the simulation domain, / is the Fourier transform of /, 
I* indicates the complex conjugate, R is the radius of the star and 2ttR^ /l"^ is the number 
of simulation domains that fit over the visible half of the stellar surface. Division by the 
Kepler finx of the simulation, F = 'Y^^w.i^i 'Y^- lij, provides for the relative power spectra of 
the flux in Eq. (Q. The we ights, Wi, and angles, /ij, of the angular quadrature are chosen 



according to the method of 



Radaij (I1880I ). which gives an optimal set of quadrature points. 



with the vertical direction included. The 0-angles that have index j, are equidistant and 
therefore have weights l/Ai'^. lijlij is the power of intensity in a particular direction, /Xj and 
(f)j. 

The power of the simulation domain is diluted by the factor 2ttR'^/P to account for 
the averaging out of uncorrelated convective fluctuations from different parts of the stellar 
disk. The underlying assumption of this formulation is that patches on the stellar surface, 
more than the distance / apart are uncorrelated. The simulations are dimensioned such 
that they each contain about a dozen major granules at any one time, and the assumption 
is therefore reasonable. 

The time series of the simulations are unfortunately not long enough to sub-divide 
the time series in order to perform ensemble averaging to reduce the noise. Instead we 
performed a running Gaussian smoothing of the power spectra and these are then fitted 
with a Harvey-like model. In the fitting process we also allow for a white noise component, 
and the two strongest p modes. Due to the rather short time series, these modes are very 
broad, and we fit them with Gaussian profiles. Only two or three modes are visible in 
each simulation, and they are constrained to have a node at the bottom of the simulation 
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box, which means they don't correspond to real eigenmodes of the star. Their frequencies, 
however, are still within the p-mode bump of real stars. 



5.2. Results on the granulation characterization 

Having fitted the granulation power spectra of the simulations of the grid to the 
generalized Harvey function of Eq. (2), we now have three granulation parameters as a 
function of T^s and log (7. For comparison with the red giants studied in this paper, we 
therefore interpolated these parameters between the simulations, to the red giants based on 
their atmospheric parameters, T^s and logg. Since the simul ation grid is irregular, we used 



linear interpolation on a Delaunay triangulation of the grid (iRenka 



19841 ). 



Top panels of Figures [TO] and [TT] show how Tes of the simulations varies with Tes and 
log g. The simulations have the same trend as the observations. We observe very tight 
correlations between T^g and log g. In the correlation between the granulation parameters 
and Teff (Figure [TO!) the stellar evolution effects are again dominant (see also Section 4.2) 
and these figures show similar structure as seen in the H-R diagrams in Figure [HI For t^s, 
the agreement between the observations and the simulations is better than a factor 2. We 
also compared Pgran from the simulations and the observations (bottom panels of Figures [10] 
and [TTj) and we note a discrepancy of an order of magnitude. 



6. Discussion and conclusions 



For ~ 1000 red giants we studied 13 months of data obtained by the Kepler mission. 
These data were processed in order to reduce the instrumental effects while keeping as much 
information as possible to preserve the granulation signal. 
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Six teams fitted the power spectra with different methods based on a Harvey model 
to estimate the granulation parameters. We first checked that the values obtained by the 
different methods were in general agreement. We notice that two branches appear for two 
of the methods that fitted only one Harvey model. In these cases, the methods are likely 
to fit for another feature present in the data. Judging from the correlation of the second 
branch with z/maxj we think this signal could also be of stellar origin. The comparison of 
the results of the red giants with the values we have for the Sun, suggests that the second 
branch would correspond to some high-frequency phenomena, such as bright points. We 
regard a further investigation of this matter beyond the scope of this paper. 

We found that the granulation time scale is proportional to while the granulation 

power is proportional to fmax^°! which agrees with theoretical scaling laws. Prom these 
results, we can use i/max as an initial guess for Tgran to fit the background. 

We also showed that stars with larger radii have larger granules to the extent that their 
surfaces are covered by a smaller total number of granules. 

Our study of the granulation power of red giants shows that larger stars present larger 
intensity fluctuations. Part of the reason for this is the smaller total number of granules 
covering the surface of larger stars, and hence less averaging out of fluctuations, compared 
to a star with many more (unresolved) granules. We also found that stars in the red clump 
have very similar values of granulation parameters. 

Pinally, we compared 3D simulations with our red giants sample and the granulation 
time scale is consistent with the observations in terms of trend with less than a factor of 
2 difference. The fact that for the observations and simulations, we have some dispersion 
in the plot Tes vs. log g can be due to several factors: 1) different stars having different 
metallicity (while our simulations have solar metallicity) , 2) the instrumental noise that 
might have not been completely removed after the data processing, and 3) the radius is not 
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necessarily determined with a sufficient accuracy for the interpolation of the simulations. 
The magnetic activity can also have an impact and could change the estimation of the 
granulation time scale by some factor, leading to different values of t^s for stars that have 
similar i/max- 
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A. Influence of the method and the initial guesses on the fitting 

We did a few tests with one of the methods used in this work (namely A2Z) to 
investigate the impact of the initial guess of the slope, a, for the following values: [1,2,3,4,5] 
on six different red giants, ranging in u^nax between 40 and 100 yuHz. To quantitatively 
compare the fits, we computed values as the mean value of the difference between the 
background fit and the heavily smoothed PDS. The best fits for the stars with i/max > 
40 /iHz were obtained for a slope of 2 and 3, with values ranging from 0.7 to 3.7, while 
for other values (1, 4, and 5), we could not reproduce the knee of the Harvey model and 
obtained only a straight line with a slope (see Figure [T2|) . For stars with z^max ~ 40 /iHz, the 
best fit is found for an initial guess of 3, 4, and 5 for the slope with values around 0.2. 

We also noticed that by fixing the value of the white noise parameter the results were 
the same except in the case of a slope lower than 3. It seems that by assuming a slope too 
large and by adding a constrain on the white noise component, the background fit becomes 
less reliable. Then we added a Gaussian to the fit to incorporate the p-mode bump. These 
values were close to the ones obtained by not taking into account the bump of the oscillation 
modes. 

We also checked how the fit converged when we decreased the initial guess of l/rgran 
to 10 /iHz. Though the code does not converge for high values of t'max, the second branch 
starts to disappear for z/max < 40/iHz. 

In the following, the fit done by A2Z with one Harvey-like model with an initial guess 
for l/xgran of 15 /iHz is chosc as the reference fit. 

We fitted two Harvey-like models using different initial values for time scale of the 
second Harvey model: 2, 5, 10 and 15 /iHz. Depending on the value of z/max, the code could 
converge or not. For all the stars, a result for the fit could be obtained when the time scale 
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value of the second model was 5 or 2 yuHz. For these stars, the between the fit and the 
PDS is smaller when the guess of the first slope is smaller than 3 and when we do not fix 
the value of the white noise parameter (around 3.6 compared to 7.5). We compared the 
dispersion of the values of t^e and of Pgran with the uncertainties for the fits that reproduced 
the knee of the granulation and found that the dispersion is related to the way we fit the 
background. Depending on the initial guess of l/rgran, we find that t^q and Pgran vary by 
~ 5% around the value found by the reference fit. The uncertainties on these values are a 
~ 6-10% around the value found by the reference fit. We conclude that the dispersion of 
these parameters with the initial guess of 1/ r^ran is of the same order of magnitude as the 
uncertainty. 

We also performed some tests on 700 red giants by fitting two Harvey-like models 
without the white noise component and with a guess for the second r of 300 /iHz. For 
these tests, a triangular smooth over 30 days was used to have enough points for the second 
Harvey like model. In Figure [131 we can notice that for these cases, Xefj is larger compared 
to the case where we use only one Harvey model. As we fit two Harvey-like models, the 
code compensates the presence of the second component by increasing the value of Tcq. 
Beside we find the same value of a in the relation TcfT=(z/max)° whether we fit one or two 
Harvey-like function while the second branch below 40 yuHz starts to disappear and to 
merge with the first one when we fit two Harvey-like functions. 

Finally, we checked that the granulation signature in the PDS was not a result of 
observing some harmonics of the p-mode bump. We removed the p modes from the PDS 
by applying a pre-whitening and fitted the background with OCT. We obtained the same 
values for the granulation characteristics confirming that there was no reflection effect. 
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Fig. 1. — Power density spectrum of a typical red giant, KIC 11618103, smoothed over 10 
bins and with the background fitting results from all the different methods. 
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Fig. 2. — Comparison of the relation between t^s and z/max for the different methods. The 
dashed line represents the fit of t^s using all the methods together. The grey symbols 
represent the second branches of OCT and A2Z. 
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Table 1. Summary of the background fitting methods (part 1). 



ID 



OCT 



SYD 



CAN 



Model 

Additional 
component 

Free parame- 
ters 

Minimization 
method 

Uncertainties 



TIT- I '^gran txt- i 



w 



Non-linear 
least squares 

Std dev. of 
parameters 



1 -I- (2TnyTiy + (27ri/Ti) 



<^i, Ti, W 



Levenberg-Markward 
least squares 

Std dev. of 
distributions from simulations 



1 + {2-KVTif 



Gaussian function 



Ti, W 



Bayesian Markov-Chain 
Monte Carlo 

Posterior probability 
of density distribution 



Table 2. Summary of the background fitting methods (part 2). 



ID 



DLB 



COR 



A2Z 



Model 

Additional 
component 

Free parame- 
ters 

Minimization 
method 

Uncertainties 



I fgran 4(7 Tgran 

1 ^ {2-KVTs,^-^f 1 + (27rl/Tgran)" 



Cr, Tgran, W 



Linear 
least squares 



C«i '^^ Tgran 



Least squares 



Std dev. of 
parameters 



. 2 

^ + 1 + (27rl/Tgran)" 



Power law 



en, Tgran, W 



Maximum likelihood 
estimator 

Inversion of 
hessian matrix 
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Fig. 3. — Granulation power obtained with the background fitting, Pgran as a function of 
^uiax- Same legend as in Fig.|2l The dashed line represents the fit of Tes using all the methods 
together. 




Fig. 4. — Granulation power Pgran as a function of the granulation time scale, r^g. The 
dashed line indicates a linear fit through the results. The second branch of OCT and A2Z 
follows the same relation. 
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Fig. 5. — Intens i ty vari ation obtained by the A2Z method as a function of log g obtained by 
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Table 3: Values of the parameter s in the power laws fitted in the Figures El [3l and HI The 
fits for OCT, A2Z, and all the methods together do not take into account the second branch. 



Method 


Teflf OC (t'max)'' 




.Pgran OC {t^sY 


OCT 


-0.82 ± 0.02 


-1.2 ± 0.02 


2.11 ± 0.03 


SYD 


-0.90 ± 0.004 


-2.12 ± 0.01 


2.37 ± 0.01 


CAN 


-0.86 ± 0.005 


-1.73 ± 0.02 


1.99 ± 0.02 


DLB 


-0.86 ± 0.01 


-2.06 ± 0.02 


2.34 ± 0.02 


COR 


-0.90 ± 0.005 


-2.15 ± 0.12 


2.34 ± 0.01 


A2Z 


-0.79 ± 0.008 


-2.09 ± 0.16 


2.39 ± 0.02 


All 


-0.89 ± 0.005 


-1.90 ± 0.01 


2.19 ± 0.01 
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Fig. 6. — Relation between the granulation time -scale, Tpff, obta i ned by the CAN method 
and stellar parameters computed as described by iKallinger et al.l (j2010al ) . 
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g. 7. — Same as Figure [6] for the granulation power, Pgran obtained by the CAN method. 



-41 - 




3.72 3.70 3.68 3.66 3.64 3.62 3.60 3.58 3.72 3.70 3.68 3.66 3.64 3.62 3.60 3.58 

log T,„ log T,„ 



Fig. 8. — For our sample of red giants, distribution of I/tcs (left panel) and Pgran (right 
panel) shown with the color code in the H-R diagram. The grey lines represent BaSTI 
evolutionary tracks computed with solar metallicity. 




Fig. 9. — Normalized emergent flux computed for a red giant simulation using the opac- 
ity distribution functions as explained in Sect. 15.11 The "lines" seen in this spectrum are 
therefore not actual spectral lines, but 1100 "giant" lines, each quantifying the distribution 
of opacity in the wavelength interval of each line (of about 20 A in the optical). The full 
calculation spans the spectrum from 860 A to 20 /im. The blue line represents the Kepler 
transmission filter. 
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Fig. 10. — Top panel: Granulation time scale, Tqs, obtained with the 3D simulations vs. Tes 
(black crosses) and interpolated to the observed sample of red giants. Results of the DLB 
method for the observations are over plotted (red diamonds). Bottom panel: Same as Top 
panel for the granulation power, Pgran- 




Fig. 11. — Same as Figure [TO] but as a function of log g. 
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Fig. 12. — PDS of KIC 11618103 smoothed over 10 bins with the comparison of the back- 
ground fitting using 1 Harvey-hke model and different values for a. The fits using a=l and 
4 are not visible because they are exactly the same as the fit using a=5. 




Fig. 13. — Comparison of the background fitting using 1 Harvey-like model (black crosses) 
and 2 Harvey-like models (blue diamonds) for a subsample of red giants. 



